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We analyze the data of purposeful rupture experiments with LOx and LH2 tanks, the 
Hydrogen-Oxygen Vertical Impact (HOVI) tests that were performed to clarify the ignition 
mechanisms, the explosive power of cryogenic H2/Ox mixtures under different conditions, 
and to elucidate the puzzling source of the initial formation of flames near the intertank 
section during the Challenger disaster. We carry out a physics-based analysis of general 
explosions scenarios for cryogenic gaseous H2/Ox mixtures and determine their realizability 
conditions, using the well-established simplified models from the detonation and deflagration 
theory. We study the features of aerosol H2/Ox mixture combustion and show, in particular, 
that aerosols intensify the deflagration flames and can induce detonation for any ignition 
mechanism. We propose a cavitation-induced mechanism of self-ignition of cryogenic H2/Ox 
mixtures that may be realized when gaseous H2 and Ox flows are mixed with a liquid Ox 
turbulent stream, as occurred in all HOVI tests. We present an overview of the HOVI tests 


to make conclusion on the risk of strong explosions in possible liquid rocket incidents and 
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provide a semi-quantitative interpretation of the HOVI data based on aerosol combustion. 


We uncover the most dangerous situations and discuss the foreseeable risks which can arise 


in space missions and lead to tragic outcomes. Our analysis relates to only unconfined 


mixtures that are likely to arise as a result of liquid propellant space vehicle incidents. 
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gas pressure Lox 
gas temperature K 
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I. Introduction 


The Challenger disaster of 1986 provoked studies of different risks that can lead to similar catastrophic events 
related to the use of H2/Ox cryogenic fuels. As was established by the Challenger investigation, the original source 
of the disaster was freezing of the O-ring in the lower section of the left solid booster and formation of a gas leak 
through the O-ring [1-3]. This leak developed into a strong jet of hot gases from the booster and caused the 
separation of the lower dome from the rest of the LH2 tank. As a result, the tank began to accelerate upwards under 
the action of the gas pressure inside the tank. The accelerating LH2 tank broke the LOx feed line in the intertank 
space and then the LOx tank’s bottom dome. The resulting LOx stream from the broken LOx feed line was injected 


into the intertank space, mixing with the GH2 jet from the rupture of the LH2 tank top dome. The mixture of fuels in 


the intertank space then self-ignited, causing disintegration of the Orbiter and the tragic loss of life (Fig. 1). 
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Fig. 1 Initiation of the first flash near the Challenger’s Orbiter/External tank forward attachment. 


The Challenger disaster represents only one possible scenario of a sequence of catastrophic events involving 
potentially explosive cryogenic fuels such as LH2 and LOx. Another such scenario has to do with an uncontained 
failure of the first stage of an LH2/LOx-liquid rocket shortly after launch. As a consequence of such an event, the 
fully loaded tanks of the second and third stages would come crashing down to the ground, violently releasing their 
entire content into the air. To predict the power of the ensuing explosion, a number of factors determined by the 
physical processes leading to the explosion have to be taken into consideration. In a situation in which the LH2 tank 
hits the ground first, the following sequence of events will take place. First, as the LH2 storage tank disintegrates 
upon impact, LH2 is ejected from the rocket onto the ground. Second, the resulting splash of rapidly evaporating 
LH2 produces GH2 and a spray of LH2 droplets in the air that are expanding from the impact location. Third, after 


some delay the rapture of the LOx tank leads to ejection of LOx and the formation of a LOx spray into the GH2-rich 


area. Direct contact between LH2 and LOx is known to lead to self-ignition of the hydrogen/oxygen gas mixture [4]. 
The energy released by the LH2/LOx combustion then vaporizes the liquid propellants, increasing their mass in the 
gas form and makes them available for further reaction, ultimately leading to an explosion. 

We note that the analysis of the underlying physical processes and catastrophic risks associated with the use of 
cryogenic fuels in rocket engines presents a number of puzzles. Up until now, the most enigmatic event in the 
sequence leading to the Orbiter’s disintegration in the Challenger disaster has been the initial formation of flames 
near the intertank section, not near the engine nozzles [1-3]. In other words, the mixture of cryogenic GH2 and 
LOx/GOx self-ignited near the intertank region. This is quite surprising, considering the fact that in the tanks LOx is 
stored at a very low temperature of about 90K and LH? is stored at an even lower temperature of only about 20K, 
while the minimum temperature of GH2/GOx mixture self-ignition at atmospheric pressure is about 850K [5]. 
Similarly, the power of the explosion following an uncontained first stage failure should depend on the degree of H2 
and Ox mixing before the ignition occurs. Clearly, only a (possibly small) part of the propellants in which they are 
well-mixed can participate in the chemical reaction. At the same time, the mixed fraction of the propellants 
sensitively depends on the time delay between the propellant release and the moment of self-ignition. On the other 
hand, the character of the explosion (a strong blast or a weak deflagration flame) should strongly depend on the 
initial density and temperature of the mixed propellants. 

To clarify the important questions about the mechanisms governing the outcome of explosions involving 
cryogens, a set of experiments with LOx and LH? tanks, the Hydrogen-Oxygen Vertical Impact (HOVD) tests were 
carried out in NASA Johnston Space Center White Sands Test Facility. We appreciate Dr. Frank Benz for providing 
us the HOVI test data. The configuration of the tanks in the tests was similar to the one in a launch vehicle, as well 
as in the External Tank of the Space Shuttle (Fig.2). 

In all HOVI tests, the LOx tank was placed above the LH2 tank. The tanks were made from an aluminum alloy 
and were insulated with 1+2 inch-thick polyurethane foam. Either the LOx tank alone or both the LH2 and the LOx 
tank were fixed on a 76 m (250 ft)-high drop tower. Then the tanks were dropped to the ground. The results of these 
tests showed some surprising characteristics of the explosions involving these cryogenic substances (see Section I 
for details). 


The purpose of this paper is to reconstruct the general picture of different types of cryogenic H2/Ox mixture 
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Fig. 2 Two main types of LH2 and LOx tanks used in the HOVI tests. The first type: two rupture devices are 
located under the bottoms of both tanks (a) and (b); the second type: one rupture device is between the LH2 
and the LOx tank (c). 


explosions based on well-know results of combustion theory [5-8] and to present an overview of the HOVI tests to 
make conclusions about the risk of a strong explosion in possible liquid rocket incidents. We provide a semi- 
quantitative interpretation of the HOVI data and analyze the ignition conditions and the parameters of different types 
of cryogenic H2/O2/N2 mixture explosions including the ones involved in the HOVI tests. We carry out a physics- 
based analysis of the general explosion scenarios and determine their realizability conditions, using the well- 
established simplified models from the detonation and deflagration theory. We uncover the most dangerous 
situations and discuss the foreseeable risks which can arise in space missions and lead to tragic outcomes. 

The purpose of this paper is to reconstruct the general picture of different types of cryogenic H2/Ox mixture 
explosions based on well-know results of combustion theory [5-8] and to present an overview of the HOVI tests to 
make conclusions about the risk of a strong explosion in possible liquid rocket incidents. We provide a semi- 
quantitative interpretation of the HOVI data and analyze the ignition conditions and the parameters of different types 
of cryogenic H2/O2/N2 mixture explosions including the ones involved in the HOVI tests. We carry out a physics- 
based analysis of the general explosion scenarios and determine their realizability conditions, using the well- 
established simplified models from the detonation and deflagration theory. We uncover the most dangerous 


situations and discuss the foreseeable risks which can arise in space missions and lead to tragic outcomes. 


We note that our analysis is limited to unconfined mixtures that are likely to arise as a result of liquid propellant 
space vehicle incidents. An independent problem is the H2/O2 mixture explosion in confined spaces, e.g. in the 
propellant feed lines. In such cases, detonation may arise as a result of the interaction of a deflagration flame with an 
external shock wave or localized obstacles [6, 9]. These effects are not analyzed in the present paper. 

Our paper is organized as follows. The main data for the HOVI tests are analyzed in Section II. The theoretic 
basis, the conditions and the parameters of detonation, deflagration, and aerosol combustion in cryogenic H2/Ox/N2 
mixtures, and also the interpretation of the HOVI test data are presented in Section III. The mechanism of 
cavitation-induced ignition of these mixtures is considered in Section IV. The conditions, parameters, and risk of the 


strong blast onset in cryogenic H2/Ox/N2 mixtures are discussed in Section V. 


II. HOVI Test Data: Conditions and Typical Scenarios of Explosions 


Most HOVI tests used two tanks, one LH2 tank and one LOx tank (Fig.2a and 2b), except for HOVI 9 and 10 
tests, in which the side-by-side double tank configuration consisting of four tanks were used (Fig.2c). The LOx and 
LH2 tanks in most tests, including HOVI 9, 13, and 14, were fixed on a 76 m (250 ft)-high drop tower. Then both 
tanks were dropped to the ground. In HOVI 2 and 5 only the LOx tank was dropped on the LH2 tank situated on the 
ground. The impact velocity was within 30+35m/sec. All HOVI tests can be divided into two groups based on the 
location of the rupture devices. The group 1 consists of the HOVI tests in which two rupture devices located under 
the bottoms of the LOx and the LH2 tank were used (see HOVI 13 and 9 in Fig.2a and b). Group 2 consists of the 
tests in which only one rupture device located between the LOx and the LH2 tank was used. HOVI 2 and 5 (Fig.2c) 
belong to this group. 

The tanks used had three typical sizes. The LH2 (LOx) tanks in HOVI 13 and 14 have diameter D=0.94 m, 
height H,= 3.84 (1.42) m, the initial fuel mass of about m2=129 (840) kg, the total volume V,=2.396 (0.817) m? and 
the gas (ullage) volume V,9 =0.553 (0.081) m’. The pressures in the tanks are pz = 1.43 atm and po,=3.15 atm, 
respectively. HOVI 9 and 10 have double tanks, with each of the LH2 (LOx) tanks having D=0.46 m, H,= 1.78 
(0.71) m, mz2=33 (176) kg, V,= 0.273 (0.095) m’, V,.o= 0.037 (0.018) m, PxH2=1.43 atm (po,=3.42 atm), respectively. 
The LOx and LH2 tanks in HOVI 2 and 5 have D=0.58 m, H,= 2.24 (0.86) m, my2=37 (189) kg, V,= 0.545 (0.185) 


m, V.o= 0.016 (0.019) m’, PH2=1.43 atm (po,=3.38 atm in HOVI 2 and po,=5 atm in HOVI 5), respectively. 


Many pressure sensors and also three film and video cameras were used to detect the explosion parameters. The 
main pressure sensors were located along three legs at ground level, 10 sensors in each (Fig.3). The main purpose of 
these tests was to obtain explosion data that would be more typical or more representative of a launch vehicle failure 


than the distributive mixture tests. 
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Fig. 3 The HOVI test site and the location of the main pressure sensors. 


Both the HOVI and the LH2/LO2 pan (dewar) tests demonstrated that the ignition always occurred and was not 
due to external sources. The HOVI test data also showed that a liquid hydrogen spill alone is not likely to self-ignite, 
because in every HOVI test with a ground cloud of hydrogen, caused by a rupture in the bottom of the hydrogen 
tank, the ground cloud did not ignite until liquid oxygen was released. The HOVI data verified the tendency for self- 
ignition of H2/Ox mixtures, because each HOVI test ignited without external assistance. HOVI test data showed that 
self-ignition occurs when GH2, GOx, and LOx mixture is available. 

In each of the test in the first group (HOVI 9, 13, 14) the LOx and the LH2 tanks were raised together. After the 
impact with the ground, first the lower rupture device (Fig. 2) broke the bottom dome of the LH2 tank, then the 
upper rupture device broke the LOx tank’s bottom dome with some delay. As a result, H2 and Ox liquid streams 
were ejected from the tanks. These streams were fragmented and the liquid droplets were partially evaporated. The 


hydrogen aerosol cloud arose near the LH2 tank bottom area (Fig.4, left). 
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Fig. 4 Dynamics of the formation of H2 and Ox aerosol clouds and the explosion for HOVI tests of the first 
(left) and second (right) types. 


The oxygen aerosol cloud appeared in the region between the tanks with the delay time tgeiq,~20-30msec. Then 
these aerosol clouds mixed partly and the explosion was observed with the total delay ftgeiay~/00msec that coincides 
with the time of the fall of LOx pieces onto the ground. The size of the observed clouds was about 2x1.4x3 m’ when 
the explosions arose. These explosions are characterized by the sensor data presented in Fig. 5. These data are 


discussed later on. 
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Fig. 5 Typical data of pressure sensors located near the explosion center for the HOVI tests of the first type. 
Data are provided by NASA Johnston Space Center White Sand Test Facility. 


In the tests of the second group (HOVI 2 and 5) the LOx tank was dropped onto the LH2 tank that was 


positioned at ground level. After the impact of the LOx tank with rupture device located between the tanks (Fig. 2c) 


the rupture device broke the top dome of the LH2 tank and, with some delay, the bottom dome of the LOx tank. As a 
result, GH2 and LOx streams escaped from the tanks into the region between the tanks (Fig. 4, right). 

The LOx stream was fragmented and the liquid droplets were evaporated partly. The oxygen aerosol cloud 
appeared with the delay time tgeiay~20-30msec. The gaseous hydrogen cloud mixed partly with the LOx aerosol 
cloud and the explosion was observed with the total delay tgejay~60ms. These explosions are characterized by the 


sensor data presented in Fig. 6. 
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Fig. 6 Typical data from the pressure sensors located near the explosion center for the HOVI tests of the 
second type. Data are provided by NASA Johnston Space Center White Sand Test Facility. 


III. Explosions of Cryogenic H2/O2/N2 Mixtures 


To analyze the experimental data and to estimate the risks of a strong explosion, we first consider detonation and 
deflagration characteristics of cryogenic GH2/GOx/GN2 mixture explosions in unconfined areas as functions of 
their temperature and composition, including the conditions for initiation of these explosions. Then, we compare 


these detonation and deflagration characteristics with the HOVI test data. 


A. Stationary Detonation Wave in GH2/GO2/GN2 Mixtures 

Detonation is supersonic combustion induced by a strong shock wave propagating directly ahead of the 
combustion front. The stationary detonation wave in unconfined reactive gas mixtures is described by the Chapman— 
Jouguet theory extended by Zeldovich, Von Neumann and Doering (ZND) [6-9]. The ZND theory is based on the 
equations of mass, momentum and energy conservation, the ideal gas equation of state, and the equations describing 
the kinetics of chemical chain reactions. The ZND model with full chemical kinetics including 21 chemical reactions 
(CANTERA [6]) was used to calculate the characteristics of the detonation waves in different GH2/GO2/GN2 


mixtures. 


The structure of the detonation wave in the GH2/GO2 stoichiometric mixture and the ZND method of its 


construction are shown in Fig. 7, where the Hugoniot curve is given by 


y (pp) p-pll 1). (1) 
y-1p p,) 2 Caria. 


and the Rayleigh line, which is tangent to the Hugoniot curve at the Chapman-Jouguet point. Here T, p, p are the 


temperature, pressure and mass density of the combustion products after the detonation wave front, and the subscript 


+ + + + + + 
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Fig. 7 A stationary detonation wave in the H2/O2 stoichiometric mixture (2:1) (ZND theory). 


T, 


mix 


“mix” denotes an initial state of the mixture before the wave front; Q = DS Y,h, (T)- 3 Y 7 ht ) is the 


reaction enthalpy, h; and Y; are the enthalpy and the mole fraction of the i-th component of the mixture. The major 


species profiles for CJ detonation in an initially stoichiometric mixture of hydrogen and oxygen (2:1) are presented 


in Fig. 8 
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Fig. 8 Major species profiles for CJ detonation in the H2/O2 stoichiometric mixture (2:1). 
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The characteristics of the detonation waves for several GH2/GO2/GN2 mixtures are shown in Fig. 9. Our 
simulations showed that, as expected, the maximum temperature T,,,,. and velocity vy, of the detonation wave 
strongly depend on the composition and weakly on the initial mixture temperature T,,,;,. Conversely, the maximum 


detonation pressure p depends strongly on Ts and relatively weakly on the mixture composition (Fig. 9). The 
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Fig. 9 The parameters of stationary detonation waves in typical H2/O2 /N2 mixtures. 


maximum pressure behaves approximately as p,,,, <(1/T,,,,)", wheren ~ 1.7. The closer the mixture composition 


to the stoichiometric H2/O2 composition (2:1), the higher the pressure and temperature of the detonation wave. The 
highest pressure (~110atm) and temperature (~3900K) in the detonation wave, i.e. the strongest blast, are achieved 
in cryogenic stoichiometric H2/O2 mixtures (2:1) at temperatures T,,,;, <100K. The blast power is determined by the 


explosion impulse PygxTimp,» Where Tin,=L/Vvg, is the blast impulse duration, L is the size of the mixed clouds, i.e. the 


radius of a hemispherical area where the mixed clouds are localized. For example, p 7. =~10* Pa-sec for the 


max ~ imp 
explosion of the stoichiometric H2/O2 mixture localized in a hemisphere of radius 2.8m. We note that the detonation 
wave front is unstable with respect to the formation of an irregular cellular structure [6, 7]. However, the main 


parameters of real detonation waves are close to those given by the ZND theory [13, 14]. 
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B. Conditions and Dynamics of Detonation Formation in Cryogenic Unconfined Mixtures 
One of the possible detonation initiation mechanisms was pointed out by Zeldovich [6-9]. It can be realized only 
in sufficiently hot explosive mixtures in the presence of a small temperature gradient when the ignition phase 


velocity in the adjacent parts of the burning mixture is higher than the detonation wave velocity ee This mechanism 


and its realization in experiments are analyzed in detail in [10, 11]. At cryogenic temperatures, however, the 
Zeldovich mechanism is not realized, since the gas mixture is too cold. 

Obviously, the detonation of cryogenic GH2/GO2/GN2 mixtures can be induced by a strong external shock 
wave with the pressure p > Pinay. Such a shock wave can be generated by a local blast. Recent experiments show that 
detonation in the stoichiometric hydrogen/air mixture enclosed in a large hemispherical envelope can be initiated by 
a blast of 10 g of C-4 high explosive located in the center of the hemisphere [12]. The detonation velocity was found 
to be 1980 m/s, which is in good agreement with the results of the ZND theory for the stoichiometric mixture. 

To analyze the possible scenarios and conditions of detonation initiation, we will use a simplified combustion 
model taking into account that the GH2/GOx burning rate is mainly limited by the initiation reactions 
H2+Ox—OH+OH and H2+Ox—HOx+H, which have the lowest rate. The rate of these reactions can be written as 


19680 
T 


Ryam = ACL) CyxCo, = Cy / Ty. Where A(T) =5.5-10" exp(— m > /mol-sec] 


)+0.7477” sep [ 


14] and 7,,, is the time scale of Cy variation. Note that this approximation is close to the one-step mechanism of 
H2 Pp p 


Mitani and Williams [15]. In this case the dynamics of the GH2/GO2/GN2 mixture combustion are described by the 


continuity equations for the molar concentration of the mixture components C: 
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t ; t 
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as well as conservation of momentum and energy 
ou 
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Here C =)'C, M. and «; are the molar mass and the thermal conductivity of the i-th component. For the analysis 


1/2 
J , 0 ~2.86-10° ea y5G, where T,, is a reference 
ge a2, wT) TM 


we use the following values: R,=8.31 
temperature. We assume that a radially symmetric explosive mixture with temperature T,,,;,=/00K at pressure D=Datm 
occupies a hemisphere of radius L,,;, and a source of ignition is in its center. Our simulations show that the 
detonation arises when the ignition source generates a shock wave with pp=80atm >p, j=55atm inside a small area of 


radius r < r=lcm (Fig. 10). 
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Fig. 10 Formation of a detonation wave in the stoichiometric GH2/GOx mixture. Initial conditions: pp=latm 
and 7,=100K for r> rp=1cm and p»=80atm and T,=100K for r < rg=1cm. The detonation wave parameters: 
Pmax=100atm, p,=60atm, Tcx=3800K, ¥4,=3000m/s. 


Parameters of detonation wave obtained: p,,,,=100atm, T,,.,=3800K and velocity v =3000m/s are close to those 


of a stationary detonation wave estimated by the ZND model (Fig. 8). If the initial pressure pp=40atm < p,j=55atm 


inside a small area of radius r < rg=lcm, the shock wave dissipates and no combustion is initiated (Fig. 11). 


5 


ow 
[=] 


Pressure (atm) 
o.6hUS8 


0.15 0.2 0 0.05 0.15 0.2 


0.1 5 
Distance (m) 


5 0.1 
Distance (m) 


Fig. 11 The distributions of pressure and temperature in the stoichiometric GH2/GOx mixture. Initial 
conditions: pp=latm and 7T)=100K for r> rg=1cm and py=40atm and T)=100K for r < r9==1cm: no detonation. 
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When the initiating shock wave has high enough temperature, T>T;,,, then detonation arises even at pp=40atm < p,; 


(Fig. 12). In this case the parameters of the detonation wave are the same as those for the case presented in Fig. 5 
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Fig. 12 Formation of a detonation wave in the stoichiometric GH2/GOx mixture. Initial conditions: pp=latm 
and 7,=100K for r> ro=1cm and po=40atm and 7,=3000K for r < r7==1cm. The detonation wave parameters 
Pmax=100atm, p,=60atm, Tcx=3800K, ¥4,=3000m/s. 


We note that the detonation of the GH2/GO2/GN2 mixture at room temperature is initiated easier than in the 
cryogenic mixture, since the pressure of the initiation shock wave can be less (Fig. 13). At the same time, the 
pressure and temperature of the detonation wave in GH2/GO2/GN2 mixtures are smaller than those for the 


stoichiometric mixture. More general conditions for detonation of the GH2/GO2/GN2 mixtures are discussed in next 


sections. 
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Fig. 13 Formation of a detonation wave in the GH2/GOx/GN2 mixture (2:1:4). Initial conditions: pp=latm 
and 7,=100K for r> 7)=1cm and pp=35atm and 7T,=3000K for r < r;=1cm. The detonation wave parameters 
Pmax=80atm, pj=45atm, Tiq.=2800K, v4,=2000m/s. 
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C. Deflagration in unconfined GH2/GO2/GN2 mixtures 

Deflagration is slow subsonic combustion mediated by heat conduction: the hot burning gas heats the adjacent 
layer of the cold gas and ignites it. There are three main processes determining the deflagration flame dynamics at 
pressures close tO Pam: (1) conductive heat transfer from the flame front to the cold mixture, (ii) burning rate 
enhancement due to turbulence, and (iii) thermal expansion of hot combustion products. The burning rate due to the 


conductive heat transfer in a quiescent gas can be estimated as 


R 
Kain = (1.5+2)m/sec (6) 


where T,,, = 310° sec is the typical reaction time of the GH2/GO2/GN2 mixtures. Turbulence accelerates the 
flame front by increasing the effective flame area. According to experimental and numerical studies, the burning rate 
increases approximately as [16]: 
Vip ~ 3-6V,) = (3+5)m/sec (7) 
Moreover, the flame front velocity increases due to expansion of the hot combustion products as 
Fam 
V pont = Vrurb a = (20+150)m / sec, 


mix (8) 
pT = Pn2Qua ~ (2000 +3900)K 


flame ~ 


C er 
P products Pp, products 


for Tix =100K+300K. 


Deflagration can be initiated by a small spark or any hot object with temperature T >T;,,, (Fig. 14). 
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Fig. 14 Formation of a deflagration wave in the stoichiometric GH2/GOx mixture: the pressure and 
temperature distributions at different times as a function of distance from the center of the hemisphere 
combustion wave. Initial conditions: pp=latm and 7,)=300K for r> ro=1cm and po=latm and 7,=1800K for r < 
ro=i1cm (in red). The deflagration wave parameters: p,,,,,=latm, T,,,,=2500K, v=28m/s. 
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Recent studies of stoichiometric H,/Ox mixture explosions show that in atmosphere the flame front velocity is vy 
=20m/sec-33m/sec and the pressure is close to latm [17]. 

The simulation results of the simplified model given by Eqs. (2)-(5) agree with the results obtained from the 
analytical estimates above and experiments [17]. The pressure in the deflagration wave is very close to the 
atmospheric pressure Pgr,. The pressure length scale is much greater than that of temperature, i.e. the “temperature 
wave’ is more localized than the “pressure wave”. This is in contrast to the detonation wave, where the temperature 


and pressure waves have the same length-scale (see Section A). The pressure in the deflagration wave falls to the 


atmospheric independently of the initial pressure of the initiating shock wave, if p < p,; (Fig. 15 and Fig. 16). 
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Fig. 15 Formation of a deflagration wave in the stoichiometric GH2/GOx mixture. Initial conditions: 
Po=iatm and 7,=300K for r> 7 >=1cm and ps=35atm and T,=1700K for r < r;=1cm (in red). The deflagration 
wave parameters: Djqx=1atm, Tyq,=3000K, v=30m/s. 
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Fig. 16 Formation of a deflagration wave in the GH2/GOx/N2 mixture (2:1:4). Initial conditions: pg=latm and 
T =300K for r> rg=1cm and po=35atm and 7,=2500K for r <7rg=1cm. The deflagration wave parameters: 
maxzlatm, Tax=2200K, v=25m/s. 
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D. Differences between the HOVI test data and the parameters of the detonation and deflagration waves 
in cryogenic GH2/GO2/GN2 mixtures 

The analysis above shows that the deflagration waves propagating in cryogenic GH2/GO2/GN2 mixtures have 
pressure p ~ Pan =latm and flame velocity Vgeq~ (25+100) m/sec, while the detonation waves have p > 40atm and 
Vaw >2000m/sec. At the same time, the sensor data showed that the maximum pressure in all HOVI tests, except for 
HOVI 9, were about (3+5) atm and the front velocities were about (660+780) m/sec (Fig. 5 and Fig. 7). Thus, these 
data conflict with the predictions for the detonation and deflagration waves propagating in cryogenic 
GH2/GO2/GN2 mixtures. We point out that in reality aerosol H2/O2/N2 mixtures, not pure gaseous GH2/GO2/GN2 
mixtures, appear in the HOVI tests. The HOVI data may be explained by aerosol combustion of these mixtures, as 


we discuss below. 


E. Formation of aerosol H2/O2/N2 mixtures in HOVI tests 

Let us first discuss the first HOVI group (Fig. 2 a), in which the impact of the tanks with the rupture devices 
results in the rupture of the bottom domes in both the LH2 and the LOx tank. The turbulent LH2 jet from the breach 
impinges on the hot ground and breaks into droplets. The rupture of the bottom dome of the LO2 tank occurs with a 
delay time tyeiay ~ 20-40msec for HOVI 13 and 14. The escaping LO2 stream brakes into droplets after the impact 
with the LH2 tank top (Fig. 4). The droplets scatter from the tank surface during ~60 msec. 
Fragmentation of a liquid stream into droplets is a complex and poorly understood phenomenon. Droplet sizes may 
vary significantly. The typical droplet radius depends on the parameters of both the liquid and gas [18]. Recent 
experimental studies of liquid jets impinging on a flat smooth surface established the following empirical correlation 


for the mean droplet radius [19]: 


2 
ee Rée= diVP1 ioe div Pr, 


Hin 2 Orn 2 


T=253*10' d, Re We (gl (9) 


air 


Here tt is the dynamic viscosity, which equals to 2 , 


air 


=1.63x10° Pasec for air, H,,,, =1.32x10° Pasec for LH2, and 


Ho, = 1.9610" Pasec for LO2, v is liquid velocity, p, is the liquid density, d), is the stream (hole) diameter, and o is 


the surface tension. Assuming v=vg=30m/sec (see Section II), we find that the typical radius of the droplets in the 
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H2 cloud is about ri, 


= 8mm and in the O. cloud is about r 
.H2 2 dr, 


a 0.5 mm. We emphasize that these values are 


almost independent of the hole diameter d),. 


The typical droplets of radius rz,and mass m= 47,1, /3 bouncing off the ground move with initial velocity 


of order vg and fly through the air with temperature T ~ 300K and pressure p = latm. The droplet velocity v is 


slowed by the air drag and is governed by: 


2 

1 1 

le S46 pg y2ar? or a =i - = SRG (10) 
dt 2 Vy at GoNVs 302 C V5 


Vv 


Therefore, the velocity and the travel distance as functions of time ¢ are equal to: 


IT 


Vv Vv 


v(t)=—"—, L =r, infie) (11) 
l+t T 
Assuming the drag coefficient for a droplet to be C<0.4, we find from Eq. (11) that the travel distance for HOVI 13 
and 14 is equal to L,~ 2.7m for the typical LH2 droplets and the explosion delay time t ~90msec, and L, ~ 2.8m for 
the typical LOx droplets for the explosion delay time t ~60msec. These values agree with the observed sizes of the 
H2 and Ox aerosol clouds (Fig. 4 and Fig. 6). 

Let us now discuss the second HOVI group, in which the rupture device is between the bottom of the LOx tank 
and top of the LH2 tank (Fig. 2c and Fig.6). In this case the breach of the LH2 tank top dome results in the escape of 
gaseous H2, while the breach of the LOx tank bottom results in the escape of liquid Ox. The released LO2 flow is 
broken into droplets after the impact with the LH2 tank top. Our estimates show that the typical radius of the 


droplets in the LOx cloud is about rj,¢,=0.5 mm and the scattering distance during the explosion delay time t 


=90msec for HOVI 2 and 5 is equal to L,~ 2.7m for typical droplets. 


F. Aerosol explosion and the interpretation of the HOVI data 

The LH2 and LOx droplets are evaporated partly by the contact with the hot combustion products and air. 
Therefore, aerosol clouds containing liquid droplets and gaseous H2, Ox, and N2 are formed. During time ty 
=10Omsec, H2 and Ox clouds partly mix and the gaseous H2 and Ox mixture ignites (see Sections A and D). As we 
noted in Section C, the temperature of such a burning mixture can be about Tyam-=3000K+3900K. To understand the 


effect of the increased GH2 and GOx masses due to droplet evaporation within the aerosols, we need to pinpoint the 
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primary heat transfer mechanism from the hot combustion products to the LH2 and LOx droplets. We note that this 
mechanism should be very efficient in order to add a significant amount of extra fuel to the mixture during the short 
explosion duration of a few milliseconds in order to have a significant amplifying effect. 

A natural candidate for the primary heat transfer mechanism is heat conduction from the hot combustion 


products through the gas phase [6]. Let us estimate the fraction of the droplet mass that would evaporate by heat 


conduction during time 7, ~3msec corresponding to the typical observed duration of the pressure spike in the 


aerosol cloud. If Lp =,/K,7 / PKg =(.5 mm is the thermodiffusion length in the gas phase, where we took 


K, > Ky ref (T; 


cael Le )* =0.1W/(m-K)as the gas thermal conductivity relative to some reference value (for 


air at room temperature), P, =lkg/ mas the gas mass density, and C, ~ 10° J/ (kg: K) as the specific heat at 


3 ; “8 
constant pressure, then the heat balance Myy,9, =4A1, TK Thame / Lp, yields Mag ~4X10™ g for an LH2 


flame 
droplet of radius r,, ,,. = 8 mm, which weighs m,, = Amr} p, /3=0.15 g. This estimate shows that only a small 
fraction (a few per cent) of the droplet mass may evaporate during the short explosion time. A similar situation takes 
place for LOx droplets of radius r,, 5, = 0.5 mm . This is due to the fact that in an aerosol formed as a result of a 


splash the droplets have rather large radii, which prevents them from being evaporated efficiently by the conductive 
heat transfer mechanism. Other situation occurs for very small spray droplets [6]. 

Another possible mechanism for heat transfer from the hot combustion products to the droplets is via radiation. 
We note that in HOVI tests the explosion was accompanied by a strong flash of bright white light similar to the one 
observed in the detonation experiments on unconfined H2/air mixtures [12]. As we will show below, this 
mechanism proves to be substantially more efficient in the case of the considered cryogenic aerosols. Heated 


combustion products, mainly water, will radiate in the infrared with the maximum intensity at wavelength 


A 1 ~(1+3) fam. According to the infrared absorption data in [20], the peak in the absorption coefficient for LH2 


Trae 


occurs at wavelength Atcorp, Ly. * 2.2 dm, which lies within this spectral range of radiation. At that wavelength the 


absorption length (the inverse of the absorption coefficient) in LH2 is eae Ly. ~ 3mm. Therefore, in this spectral 
range of the infrared the droplets are opaque and are able to absorb a significant portion of the incoming radiation. 


Furthermore, since Lica LH2 18 comparable to the droplet radius, the heat from radiation will be absorbed in the 
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droplet bulk, raising the temperature of the liquid phase inside without significant evaporation at the droplet 
interface. Thus, the radiation heat may quickly raise the LH2 and LOx droplet temperature to the critical temperature 
T., resulting in an explosive vaporization of the entire droplet. This will greatly enhance the aerosol combustion. 


The evaporation time for a droplet subject to radiative heat transfer can be estimated via the following heat balance, 
C,p,(f.-1,)| = 2 |=n'oFt,.4ar 12 
iPr f c 1) Re Var ca | oO flame AVL evap 2 ( ) 


where o=5.67x10°W/m’/K* is the Stephan-Boltzmann constant and n is the absorption efficiency. We will 


estimate 7 = 0.1+0.5, assuming that a significant portion of the radiation is absorbed by the droplet. We note that 


a more precise estimation of the absorption efficiency requires a detailed analysis of a very complex problem of 
light emission by the combustion products (water) and propagation through a highly heterogeneous aerosol mixture, 
which is beyond the scope of this paper. Here we limit ourselves to semi-quantitative estimates aimed at capturing 
the orders of magnitude of the quantities of interest and allowing us to fit the experimental data. As follows from 


Eq. (12) 


T, , a = CP, ce mele (13) 
. 37 OT, 


ane arr T 
flame 


where 7;=3500K is the characteristic flame temperature. The evaporation time strongly depends on the flame 
temperature Tyq,.. Taking n 2 0.25 and the droplet radii from Section E, we have f,,,,< 10msec both for the LH2 


and the LOx droplets at T =3500K. Of course, t.,,, depends on a number of factors and may vary in the range from a 
few milliseconds to tens of milliseconds. 

Importantly, the vaporization of the LH2 and LOx droplets results in an abrupt increase of the combustible gas 
density and, therefore, leads to a strong buildup of pressure in the combustion products. The total mass of the burned 


burned 


= : * burned __ _ burned 
GH2,My) 0 = Mgyx + My7 drop > 18 controlled by the available mass of GOx, mg" = Mo, + Moy drop = 8M 2 


formed by the initial gaseous Ox and the evaporated LOx droplets, mixed with GH2. The product mass (water) is 


burned burned burned 


Mo. +M,5. =9m,,;— and the pressure of the aerosol combustion wave can be estimated as 


Pac = 9%. (Pour + Our" )+ >) Ri, Tame’ (14) 
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where ,,R, are the mass density and the gas constants for the i-th unburned gas component, 


droplet 
H2 


= mi a / Vi. , where V, =27 bi / 3 is the volume of the hemisphere in which the aerosol cloud is mixed, 
Pacw 18 the pressure in the aerosol combustion wave. Assuming 7.,,, 18 equal to the maximum pressure P,,,q, that is 
determined by the sensor data (Fig. 5 and Fig. 7) we can estimate the density of the droplets from Eq. (14) as 
Rp. 
2 is ma 8p max 
ORicy. | OReal 


flame 


droplet 


Por 


droplet a 8 Pp max 


=8 OF Dax >> Patm = atm (15) 
Pr OR, of for p Pat 


— Pou. 


ame 


To simulate the aerosol combustion we use Eqs. (2)-(5) in which we replace Eqs. (2) with Eqs. (16): 


oC Copies T 
ao V: (uC,,>) = —2Cy.Co, A(T) + ant : 


evap 1 


droplet T i 
C09 (uC y,) = —Cy CoAT) + (z ; (16) 
1 


4 4 
droplet droplet droplet droplet 
[Sp Ginn OF (F) Se glee Os (F) 
> 


dt Fite holy dt %. T, 


evap 


evap 


where C. is the molar concentration of the i-th component in the mixture. 


As we noted above, the first HOVI group (HOVI 13 and 14) is characterized by the formation of both Ox and H2 


aerosol clouds forming due to the rupture of the bottom domes in both tanks. These clouds first mix and then the 
explosion occurs. To simulate this aerosol explosion, we take into account that Cr >(Qand Cos >O and have 


the initial values giving the pressure in the aerosol combustion wave which coincides with p,,,, in the explosion 


wave measured by the sensors (Fig. 5 and Fig. 7). For example, according to Eq. (15) the sensor reading of 


Pmax~4atm in HOVI 13 corresponds to an aerosol combustion wave with C4" =20mol/m (p93? =0.04kg /m’) 


and C4" = 10mol /m?> (p%""' = 0.32kg / mr’) (see Fig. 17). 
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Fig. 17 The distributions of pressure and temperature during aerosol combustion in the mixture (2:1:4) with 
the evaporation time Cig 1ms and the droplet densities p;.=0.04 kg/m’ and po,=0.32kg/m*. Initial 


conditions: pp=latm and 7)=100K for r> rg=1cm and py=latm and T,=2500K for r < r9=1cm (in red). The 
wave parameters: Djqx=4.2atm, Ting,=3300K, v=600m/s. 


The second HOVI group (HOVI 2 and 5) is characterized by the formation of an Ox aerosol cloud and 


a gaseous GH2 cloud arising from the rupture of the LOx tank bottom dome and of the LH2 tank top 
dome. In this case Ce "' = (and the initial concentration C yin the mixture greatly exceeds C,, . 


Figure 18 shows the formation of the aerosol combustion wave with the parameters close to those of the 


HOVI5 explosion (Fig. 6). 
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Fig. 18 The distributions of pressure and temperature during aerosol combustion in the H2/O2/N2 (3/1/2) 
mixture with Le IMS, pro = 0 and po, = 0.02kg/m’. Initial conditions: pp=latm and 7,=100K for r>r,=1cem 


and a hot object with T=2500K (in red). The wave parameters: Pjgx=5.4atén, Tngx= 3200K, v= 720m/s. 


Comparison of the parameters of the simulated aerosol combustion waves with the sensor data are given in Table 1. 


This table summarizes our findings regarding the cryogen masses involved in different representative HOVI tests, as 
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well as the main characteristics of the blast. Specifically, column | lists the numbers of the typical HOVI tests for 
which the analysis was presented, column 2 shows the masses of mixed LH2 and LOx aerosols found from the 
numerical simulations of Section F to produce a blast wave with the peak pressure observed experimentally, column 
3 lists the masses of LH2/GH2 and LOx that escaped from the tanks prior to explosion computed in Section G, 
column 4 lists the peak pressures recorder by the sensors, column 5 lists the blast wave velocity measured 
experimentally, column 6 lists the values of the calculated blast wave velocity based on the shock wave speed at 
pressure shown in column 4, column 7 lists the duration of the dynamic pressure impulse measured by the sensors, 
and column 8 lists the calculated values of the blast duration computed as the ratio of the aerosol cloud radius to the 


shock wave velocity. 


Table 1 Main parameters of aerosol combustion and the sensor data 


AeGsol Measured 
mass of ees Maximum | Experimental | Calculated blast duration | Calculated 
mixed H2 and Ox | pressure blast wave wave velocity of the duration of 
H2 and mass (kg) from velocity (m/sec) blast the blast 
Ox (kg) at sensor (m/sec) 7 (1 4 (DE i wave wave 
- fie Rea tOcm | data (atm) v=AL/ At cee 2700 (msec) tT=Liv 
> 66 LH2 : oe re = ~ 
313 LOx 3.3+4.2 = 740+825 = 760 23.5 23.5 
3.2+5.4 = 780 = 785 =3.4 =3.4 
2.5+3 = 660 = 660 =4.0 =4.0 
80+110 2625+2966 2500+2928 =0.7 =0.7 


One can see from these data that only a small fraction of the escaped cryogens mass participates in the explosion 
in the form of aerosol. Also, the blast wave velocity and duration agree very well with the numbers based on the 


peak pressure and the size of the aerosol cloud. We note that deflagration wave acceleration in aerosol H2/O2/N2 


>Ims. However, for larger values of T. 


mixtures also occurs for lower droplet evaporation rates, i.e., for 7, eiap 


evap 
one obtains very similar combustion waves, provided that initial droplet mass is higher. For example, for the 2:1:4 


mixture and ‘ane =5msec one finds py =0.07kg/m* and po,=0.56kg/m*. Let us point out that an aerosol 
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deflagration wave can accelerate and transform into a strong detonation wave for large enough aerosol densities and 


sizes of the mixed H2 and O2 aerosol clouds. This effect is considered in Section G. 


G. Estimates of the escaped H2 and Ox masses before the explosion 


The aerosol combustion is determined by the H2 and Ox masses inside the overlapping area of the H2 and Ox 


aerosol __ aerosol 3 aerosol aerosol 3 : = 
aerosol clouds. These masses are m= 2. (27/ 3)L ne ape (22 /3)L,,,- The typical size of the 


overlap area is Ly, ~2m. Obviously the H2 and Ox masses escaped from the ruptured tanks during the short delay 
time tyclay<100msec before the explosion have to be greater than the aerosol masses obtained in Section F. Our 
estimates below show that this condition is fulfilled for all HOVI tests. Indeed, the escaped LH2 mass in the first 


group of the HOVI tests can be estimated from the equations of momentum, energy and mass conservation for the 


LH2 flow: 
2 2 f 
PLY _ PY Pp _| Veo c 
Be 5 Pe ee ae 
Po g C, (17) 


dh, =—dh, = = vdt, dV, =-S,dh, = S,vdt =-dV, 


0 


where S pao R, is the cross-section area of the hole in the tank bottom. It follows form Eq. (17) that 


1/2 


7 1/2 
dV, ay V0 — Patm 4 PLY gos 22) (18) 
, h 
dt V, Po 2Po L 
PY P 
For LH2 we have 9, = 70kg / m Y= CIC, = 1.4, p, = 1.45amm, p _,, = latmand the condition £0 < am 
2P) Py 
is valid for the impact velocity vo< 45m/sec. Therefore, according to Eq. (18) 
2 -lly 
* 
Vous =,( 2a _ PL i : (19) 
Py = 2Po 


It then follows that mj. = Pyi2V, max ~ 66kg for HOVI 13 (V,, = 0.553 m’) and mye" = 2kg for HOVI 14 
(V,) =0.176 m’)- The volume of the escaped LH2 equals to its maximum value V_, max for time 20ms < faeiay- Vina 


does not depend on the ignition delay time f4./,. and the cross-section S, =a, of the hole in the LH2 tank top. We 
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see that Wee = 66kg(HOVI 13) or 21kg(HOVI14) 1s much more than the necessary mass of the H2 aerosol in the 


overlapping H2 and Ox clouds: mor" = pererlery @orlet — 0.6 7kg, Vir"! = (20 / 3) L; 


H2 H2 overlap 


for par =0.04kg /m?* and L 


overlap 


= 2m. For Ox, we have 9, =1141kg /m,y= C ,/C, =1.4, 


2 
Do = 2.1atm, P im = latm , and the opposite condition PLY0_ 5, Pam ig valid for the impact velocity vg>14m/s. 
2D Po 
Analysis of Eq. (18) shows that mee > 15kg for the ignition delay time tgeigy > SOmsec and R;,> 10cm. At the 


same time mi" = phoeyroret ~ 5 Akg < mi" . Moreover, there is additional oxygen in the atmosphere. 


The second group of the HOVI tests is characterized by the rupture device located between the tanks. In this case 
the impact results in the rupture of the bottom of the LOx tank and the top of the LH2 tank. As a consequence, liquid 
Ox and gaseous H2 escape from the tanks. The escaped LOx masses are determined by Eq. (18) and their values are 


presented in Table 1. The escape dynamics of GH2 is given by [21] 


dp, 
Meo dt 
/y 
o-[2 (20) 
P 
fese max (0) - 
My =S8, ; j(t) dt, m5 a ve 
0 


The values of the escaped GH2 mass for the delay time fyejq,~ 100 msec are presented in Table 1. These parameters 


were taken into account in other calculations. 


H. Deflagration-to-detonation transition in aerosol mixtures 

As we noted above, deflagration in the explosive H2/Ox mixtures may be initiated by a small spark or any hot 
object with temperature T > Tj, (Fig. 14). Due to the presence of aerosols the deflagration wave accelerates and its 
pressure increases (Section F). When the pressure exceeds the value ~p,, a strong detonation wave arises. The 
characteristics of this wave are determined only by the parameters of the H2/Ox mixture and do not depend on the 
initial condition triggering this explosion. Thus, the greatest explosion risk is posed by the aerosol H2/Ox mixtures. 
These mixtures are easily ignited and can produce an especially strong explosion when the mixtures are close to the 
stoichiometric composition. This effect is realized at relatively high droplet densities of both H2 and Ox 


components (Fig. 19), or-even LOx alone (Fig. 20). 
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t=0.63ms t=0.7ms 


Fig. 19 Deflagration-to-detonation transition in the aerosol H2/O2 mixture (~2:1) with the evaporation time 
= lms and the droplet densities p;,.=0.08 and po,=0.64kg/m3. Initial conditions: py=latm and 7,=100K 


for r> r0=1cm and T=2500K (in red). The wave parameters: Pj, =102atm, Tngx= 3800K, v =3000m/s. Inserts 
show early stages of the processes. 


This leads us to an important conclusion: a strong detonation explosion can be initiated by a spark or a hot object 
in an aerosol mixture, not just by a strong shock wave, as is the case for cryogenic gaseous mixtures (Section B). We 
note that the maximum pressure of the aerosol detonation wave can exceed that of the detonation wave in gas 
mixtures (compare Fig. 11 with Fig. 19). The aerosol detonation scenario is apparently realized in the HOVI 9 test, 


which resulted in an unusually strong blast (see detail in Section E). 
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Fig. 20 Deflagration-to-detonation transition (3:1:0) in the aerosol H2/O2/N2 mixture with 7, =Ims, 


evap 
Pr2=0 and po,=0.48k¢g/m3. Initial conditions: pp=latm and T)=100K for r> r0=1cm and T=2500K (in red). 
The wave parameters: Pq, = 100atm, Ting, =3800K, v =3000m/s. Inserts show the beginning stage of the 
processes. 


Let us point out that the effect of deflagration-to-detonation transition in aerosol H2/O2/N2 mixtures discussed 


above also occurs at lower rates of droplet evaporation, e.g., when the evaporation time Tovap > IMs - However, for 
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larger values of T,,, ap higher droplet mass is necessary. For example, for stoichiometric aerosol H2/O2 mixtures one 


needs Pr2=0.32kg/m* at Tovap = 5ms. 


IV. Cavitation-induced ignition of H2/O2/N2 mixtures 


One of the puzzling questions in the studies of cryogenic explosions is the mechanism of ignition in the 
cryogenic H2/O2/N2 mixtures [22]. Multiple tests, including the HOVI tests, give clear evidence that the ignition in 
these mixtures is not due to external sources. GH2 released alone does not explode, despite of the large amount of 
oxygen present in the air. Ignition occurs when gaseous or aerosol cryogenic hydrogen and oxygen, as well as liquid 
oxygen (LOx) are available. The tests showed that self-ignition of the cryogenic mixtures is realized always for the 
short time less than 0.1s after GH2 and GOx flows are mixed with a turbulent LOx stream. This condition is 
fulfilled in all HOVI tests. Below we discuss a possible mechanism of self-ignition of cryogenic H2/Ox mixtures 


that can be realized when GH2, GOx, and LOx flows mix. 


I. Cavitation-induced Detonation of H2/O2/N2 Mixtures 

We propose a cavitation-induced mechanism which may explain self-ignition in cryogenic H2/Ox mixtures. 
Cavitation refers to the formation and compression of vapor bubbles in a liquid under the action of the applied 
pressure difference. Due to the inertial motion of the liquid, this process results in a rapidly collapsing bubble and an 
increase in the gas temperature and pressure inside the bubble, producing a strong shock wave [23, 24]. In the 


simplest case, the dynamics of the bubble radius is described by the Rayleigh-Plesset equation [22] 


R, 


2 
aR, 5) At, AR, | 20 _ Dy (Ry) Pin (21) 


ai? ° 2. dt p,R, dt p,R Pp, 
where p,(R,) and p,,, are the pressures inside the bubble and in the liquid far from it, respectively. A pressure 


jump p,(R,)— p, between the LOx and the bubble can arise from the turbulent mix of gaseous H2 and Ox and Ox 


liquid stream or due to relatively weak initiating shock wave arising from in the LOx piece containing the bubble as 
a result of impact of this piece on a solid object. Cavitation-induced ignition of H2/Ox mixtures may be caused by 
the second strong shock wave generated by a collapsing bubble near the LOx surface which propagates across the 


GH2/GOx-LOx interface and ignites the gas mixture (Fig. 21). 


27 


GOx/GH2 Hot flow injection 


Rate of 
liquid 


Collapsing 


GOx/GH, bubble 


Fig. 21 Cavitation-induced ignition: collapse of a vapor Ox bubble in liquid Ox under the action of a pressure 
differential p;-p, (a); the bubble can contain only Ox vapor or a GOx/GH2/GN2 mixture (b); generation of 
the shock wave and injection of hot gas from a bubble collapsing near the LOx - GOx/GH2 interface (c). 


To estimate the gas parameters in the collapsing bubble, we will assume that the bubble contains Ox vapor and 


maybe a small portion of GH2. The analysis carried out in [23, 24] show that we can assume that the Ox vapor 
pressure due to active condensation iS closed to the pressure of saturated Ox vapor, and the compression of GH2 in 


the collapsing bubble is an adiabatic process. In other words, we can write the pressure in the bubble as [23] 
37 
P(t) = Ps (Ty) + Peo (Ryo / R) (22) 


Here R,,. and p go are the initial bubble radius and the initial partial pressure of the H2 gas in the bubble, p, (7; ) is 


the saturation pressure of the Ox vapor at temperature T equal to the temperature T, at the liquid-vapor interface. The 


time of the bubble collapse can be estimated by [22] 


R2 1/2 

te = 0.915 Pi. (23) 
Pi. i Ps; 

The results of the numerical analysis of Eqs. (21)-(22) are presented in Fig. 22. These results show that the collapse 

time is close to that given by Eq. (23). We also see that under the action of the pressure difference exceeding 0.2 atm 

the initial bubble radius Rp=2mm goes to its minimum value R,,,;,~0.1mm, and the temperature and pressure in the 


collapsing bubble can exceed 500 atm and 2500K, respectively. Our simulations based on Eqs. (2)-(5) show that a 


localized shock wave of lower intensity is sufficient to induce detonation both in stoichiometric H2/Ox (Fig. 23) and 
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Fig. 22 Collapse of a GOx/GH2 bubble of initial radius R,,=2mm in LOx. Initial conditions: 1 - p,- 
Po2=0.75atm, po .=0.25atm, p;=0.015atm; 2 - p,-po2=0.5atm, po,=0.5atm, p7=0.0085atm; 3 - p;-po2=0.25atm, 
Po2=0.75atm, p;=0.003atm. The parameters are: 1 - Dyg=510atm, T,=1450K, Rpnin=0.165mm, collapse time 
teottapse=9.23MS: 2 = Prax =460atm, Tyyq.=1850K, Rjnin=0.15mmM, feonapse= 0.284mS; 3 - Pngx = 450atm, Ronin= 
0.1mm, T,,¢:=2500K, teorapse= 0.4ms. 


H2/Ox/N2 mixtures (Fig. 24). We note that in real situation the super-hot and super-compressed O, H, OH species 
are formed in process of mixture burning [6] in the bubble (see Fig. 8). These species will be ejected from the 


bubble into the space above the LOx surface and easily ignite the GH2/GOx mixture situated under it. 
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Fig. 23 Cavitation-induced ignition of a stoichiometric H2/Ox mixture (2:1) with the mixture temperature 
T=100K and pressure p=1latm for r>0.15msec. Initial cavitation condition: temperature 7;=1500K and 
pressure po=350atm for r<0.15mm (in red). The detonation wave parameters: p,,,,,=100atm, p,=60atm, 

T nax=3800K, ¥4,=3000m/s. 


The main requirement for the cavitation onset is a fast jump of pressure between the liquid and the vapor bubble. 


This jump may be of different nature. We list here possible scenarios of cavitation-induced ignition in H2/O2/N2 


mixtures: 
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Fig. 24 Cavitation-induced ignition of a H2/Ox/N2 mixture (2:1:4) with mixture temperature T=100K and 
pressure p=latm for r>0.15mm. Initial cavitation condition: temperature 7,=1500K and pressure py=350atm 
for r<0.15mm (in red). The detonation wave parameters: p,,q,=82atm, p,;=45atm, Tyq,=2800K, v4,=2000m/s. 


1. Collapse of a_ rarefied vapor bubble. Rarefied vapor bubbles can arise in the liquid as a result 


of turbulent mixing of GOx and LOx streams and collapse in it. 


2. Collapse of a vapor Ox bubble with admixed GH2. In this case ignition can be intensified due to the possibility of 
a local explosion inside the collapsing bubble because of the high temperature and pressure in it. Along with the 
generation of a strong shock wave, the super-hot and super-compressed O, H, OH species can form in the bubble 


and be ejected from it into the space above the LOx interface and easily ignite the GH2/GOx mixture. 


3. Formation of bubbles with chilled surface. The surfaces of large LOx pieces can be cooled by very cold GH2 
flows. The Ox vapor bubbles with a thin surface layer of cooled liquid and, probably, with admixed gaseous H2 as 
well, can rapidly form as the result an of impact of two LOx pieces. Due to the low surface temperature the pressure 


in the bubble quickly drops because of intense vapor condensation inducing bubble collapse. 


4. Injection of an LH2 droplet into LOx. The impact of the LH2 and LOx jets with the ground results in turbulence 
and fragmentation of the liquids into droplets. A cold LH2 droplet (with T= 20K) can be captured by a “hot” LOx 
piece (T=90K). The LH2 droplet will then evaporate very quickly and the pressure inside the bubble will quickly 
grow and can become much greater than the pressure p; in LOx. As a result, the bubble radius will increase. Due to 
the inertial motion of the liquid the bubble expansion can lead to the final pressure much less than p;. Afterwards, 


the bubble starts collapsing and the pressure and temperature of the GH2/GOx mixture inside the bubble can become 


very high and initiate a strong shock wave. 
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These scenarios were confirmed by our preliminary numerical simulations. A detailed analysis of these results is 


beyond the scope of this paper [25]. 


J. Cavitation-induced deflagration of inhomogeneous H2/O2/N2 mixtures 

H2/Ox/N2 explosive mixtures containing gaseous and/or aerosol hydrogen, oxygen, and nitrogen and LOx may 
have strongly inhomogeneous. The collapse of a bubble near the LOx interface can generate a strong shock wave of 
high pressure and temperature, inducing a detonation wave in a hydrogen-rich area located near the LOx surface. If 
this area is surrounded by a region with low H2 density (Fig. 25), then the detonation wave will rapidly dissipate in 


this region and the temperature wave can initiate a deflagration in the mixture at the other side of the hydrogen- 


depleted region. 
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Fig. 25 Transition from local detonation to deflagration: a sketch of the propagating blast and heat waves in a 
GH2/GO2/GN2 mixture generated by a collapsing bubble. 


Our simulations confirm this conclusion (Fig. 26). We see that the initial parameters of the explosion wave change 
from pressure p=350 atm to 1 atm, from temperature T=3000K to 2600K, and the front velocity from v=1200m/sec 


to 30m/sec in the forming deflagration wave. 
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Fig. 26 Transition from a cavitations-induced local detonation to deflagration of a GH2/GOx mixture with 
T=100K and p=1 atm. Initial condition: T=3000K, p=350atm inside area r<r,;=0.1mm. The parameters of the 
deflagration wave: T,,,,,=2600K, p=1 atm, v=30m/s. 


V. Conclusion: conditions and risks of the onset of strong blasts in the cryogenic mixtures 

We believe that proposed cavitation-induced ignition mechanism and considered scenarios of its realization of 
are important for understanding conditions and risk of explosion of cryogenic H2/Ox fuel used in liquid rockets and 
other vehicles. 

The main results of our analysis are summarized in Fig. 27 describing possible scenarios and conditions of 
different combustion types for cryogenic H2/Ox/N2 mixtures. The combustion type is determined by the 
composition of the mixtures, the degree of their mixing, and also by the ignition mechanism. The strong explosion 
(detonation) can arise in gaseous GH2/GOx/GN2 mixtures when they are well-mixed and the ignition is induced by 
a strong shock. The closer the mixture to the stoichiometric composition, the stronger the detonation. Any other 
ignition source, e.g. by a spark or a hot object, results in slow combustion (deflagration) in gaseous GH2/GOx/GN2 


mixtures. 
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Fig. 27 Diagram of possible cryogenic H2/Ox explosion scenarios. 


Self-ignition is realized when GH2/GOx/GN2 flows mix with a turbulent LOx stream. Any mechanism of 
ignition of partially mixed aerosol H2/Ox mixtures with relatively low droplet concentration leads to aerosol 
combustion, accelerated deflagration, and is characterized by the overpressure of several bars. Well-mixed H2/Ox 
aerosol mixtures with high droplet concentrations are most dangerous: any ignition mechanism, including self- 
ignition, in contact with relatively large LOx pieces results in a strong explosion. The maximum pressure in such an 
explosion may exceed 100 atm for stoichiometric mixtures. 

The explosive power depends on the pressure and the size of the area of mixed H2 and Ox aerosol clouds. A 
strong detonation arises when the mixture is well-mixed inside a large enough volume. The blast intensity may be 
characterized by the dynamic pressure impulse PyaxTimp, Where Timyp=L/vay is the blast impulse duration and L is the 
size of the mixed clouds. This quantity may be used to determine the TNT equivalent of the explosion [4]. The value 
of L depends on the mixing time of the aerosol clouds, i.e. the time of delay t4.j,, between the tank breach and the 
explosion. The relatively small value of tyejqy~ (30+60)msec in the HOVI tests corresponds to the case when the 
breaches occur in the intertank area and gaseous H2 and liquid Ox are injected into the intertank area from the upper 


ruptured LH2 tank dome and lower ruptured LOx tank dome. Such a situation was realized by accident in the 
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Challenger disaster [1-3] and purposely in the second group of the HOVI tests (Fig. 2c). The escaped GH2 and LOx 
flows very quickly (but poorly) mixed together and self-ignited. Similar events can occur in the case of the bulkhead 
failure [26]. 

The aerosol combustion, accelerating the deflagration combustion, is likely to arise in the cases with the 
maximum blast pressure from one to several bars. Here the risk of a strong explosion is intermediate. Different 
situation may occur when the bottom domes of both tanks are ruptured, as can be seen from the first group of the 
HOVI tests (Fig. 2a and 2b). Here the delay can be longer due to a much larger area separating the H2 and Ox 
aerosol clouds and the large LOx pieces. Such an event was probably realized in the HOVI 9 test. This test used 
doubled tanks (Fig. 2b) with the rupture devices placed underneath each tank. Following the rupture of the LH2 
tank, its thermal insulation detached and shielded the H2 aerosol cloud from direct contact with the turbulent LOx 
stream. This resulted in a long delay time tgeiay~ 200 msec between the rupture and the cavitation-induced ignition. 
During this time, the expanding H2 and Ox aerosol clouds had a chance to form and mix together. As a 
consequence, a nearly stoichiometric H2/Ox mixture confined to a hemi-sphere of radius of about 2m was formed. 
The self-ignition of this mixture induced the strongest explosion (aerosol detonation), characterized by the 
maximum pressure exceeding 100 atm and the duration of about Imsec (Fig. 24). We attribute this strong explosion 
to a chance occurrence of the complex breaking dynamics of the insulating foam detaching from the double LH2 
tank. Such a strong explosion was not observed in all other 13 HOVI tests, including the HOVI 10 test, which had 
the same double tanks as HOVI 9. In the case of HOVI 10 the value of tg.j4, was very small and pjq,=1.36 atm. 

An event in which the contents of the LH2 and LOx tanks is scattered over a relatively large area is probably the 
most dangerous from the point of view of the risk of the strongest explosion. Indeed, in this case the escaped H2 and 
Ox liquids will have large enough time to evaporate and to generate GH2 and GOx aerosol clouds that will have 
time to mix together into a large enough area before self-igniting upon contact with the ejected turbulent LOx 
stream. We emphasize once more that self-ignition of cryogenic H2/O2 mixtures is always realized when GH2 and 
GOx flows are mixed with a turbulent LOx stream. 

We acknowledge Dr. Frank Benz from NASA Johnston Space Center White Sands Test Facility for providing us 
the test data for Hydrogen Oxygen vertical impact tests performed at the facility. We also acknowledge Dr. Benz for 


interesting discussions. 
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